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Abstract 

Solitonic scalar field configurations are studied in a theory coupled to gravity. It is found 
that non-topological solitons, Q-balls, are present in the theory. Properties of gravitationally self 
coupled Q-balls are studied by analytical and numerical means. Analytical arguments show that, 
unlike in the typical flat space scenario, the size of Q-balls is ultimately limited by gravitational 
effects. Even though the largest Q-balls are very dense, their radii are still much larger than the 
corresponding Schwarzschild radii. Gravity can also act as a stabilising mechanism for otherwise 
energetically unstable Q-balls. 
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1 Introduction 



Various field theories can support stable non-top ological solitons p], Q-balls [^], which are co- 
herent scalar condensates that carry a conserved charge, typically a U(l) charge. Due to charge 
conservation, the Q-ball configuration is the ground state in the sector of fixed charge. Q-balls 
may have physical importance because the supersymmetric extensions of the Standard Model 
have scalar potentials that are suitable for Q-balls to exist in the theory. In particular, lepton 
or baryon number carrying Q-balls are present in the Minimal Supersymmetric Standard Model 
(MSSM) due to the existence of flat directions in the scalar sector of the theory [||, |], H, ||. 

Q-balls can be cosmologically significant in various ways. Stable (or long living) Q-balls are 
natural candidates for dark matter H and their decay offers a way to understand the baryon 
to dark matter ratio as well as the baryon asymmetry of the universe 0, ||. Q-balls can 
also protect the baryons from electroweak sphalerons Q and may be an important factor in 
considering the stability of neutron stars 0]. 

For Q-balls to be cosmologically significant, one needs to have a mechanism that creates 
them in the early stages of the evolution of the universe. Q-balls can be efficiently created in 
the early universe from an Affleck-Dine (AD) condensate [[|, [5]. This process has been studied 
by numerical simulations §§-|ll] where both the gauge- and gravity-mediated SUSY breaking 
scenarios were considered. In both cases, Q-balls with various charges have been seen to form 
from the condensate. Also collisions of Q-balls have been considered in various potentials [l2| ], 

A large Q-ball is a very dense object and hence its coupling to the gravitational field may be 
remarkable. In the present paper we study these aspects of Q-ball solutions i.e. we are interested 
in the gravitational coupling of scalar configurations. Some such considerations have also been 
studied with reference to Q-stars [|l3|]. 

The present paper is organised as follows: In Chapter 2 we briefly review the properties of 
Q-balls in flat space and study analytically the gravitationally coupled system. In Chapter 3 
special attention is given to gravitationally coupled thin-walled Q-ball configurations. Numerical 
results are presented in Chapter 4. The paper is concluded in Chapter 5. 



2 Theoretical considerations 
2.1 Preliminaries: Q-balls in flat space 

Consider a field theory with a U(l) symmetric scalar potential, V(|0|), with a global minimum 
at <p = 0. The complex scalar field eft carries a unit quantum number with respect to the U(l)- 
symmetry. The charge and energy of a field configuration ^ in D dimensions are Jl| 

Q = - f (<f)*<j)-W*)d D x (1) 



i 

and 



E = JiW + m' + V^W^x. (2) 

The Q-ball solution is the minimum energy configuration in the sector of fixed charge. A Q-ball 
is stable against radiative decays into 0-scalars if condition 

E < mQ, (3) 

where m is the mass of the (/>-scalar, holds. It is then energetically favourable to store charge in 
a Q-ball rather than in the form of free scalars. 

Finding the minimum energy is straightforward and the Q-ball solution can be shown to be 
of the form |, g 

0OM) = e*"V(r), (4) 



where f>(x) is now time independent and real, u is the Q-ball frequency (for most cases \ui\ £ 
[0,m]). The charge of a Q-ball with spherical symmetry in D-dimensions is given by 



Q = 2oj J <t){rfd D r (5) 
and the equation of motion at a fixed uj is 

dr 2 r dr d(j) 2 

To find the Q-ball solution we must solve (|6|) with the boundary conditions 4>'(0) = 0, </>(oo) = 0. 

The solution of Eq. (g) depends crucially on the shape of the potential. Flat potentials, 
which are essentially constant over wide range of field ip give rise to the relation E oc Q 3 / 4 Ji~4] ], 
whereas in a wide class of more general potentials the energy-charge relations is linear, E oc Q 
|l|, E|. However, any Q-ball solution in a potential which is bounded from below, asymptotically 
increases faster than ip 2 and has a global minimum at the origin, approaches the thin-wall limit 
when uj becomes small enough. At the thin-wall limit any potential then leads to the linear 
relation between energy and charge. 



2.2 Field equations for gravitationally coupled Q-ball 

Supposing that the space time has a maximally symmetric subspace S 2 , i.e. the configuration is 
spherically symmetric, the metric tensor can be written in the form 



g^u = diag( J B, -A, 



-r 2 sin 2 6>), 



(7) 



where A and B are positive functions of time and radial coordinate [15|. A complex scalar field 
(f> coupled to gravity is then described by the action 



S[(/>, A, B] = l — I d 4 xVABTZ+ I d 4 xVAB ( d, 

16-7TG J J v 



V 



(8) 



where 1Z is the Ricci curvature scalar, V(4>^(j)) is the scalar potential and G = 1/Mpi- By varying 
the action one obtains the well known Einstein equations together with the equation of motion 
for the scalar field. They read for the spherically symmetric system in the contracted form as 



B" B' (A' 
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B^_ 
2B 
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AA\A 




(9) 



where the energy density p and pressure P are given by 

L>\2 



p - 
p 



k'\2 



V. 



(10) 



The total energy E'tot of the Q-ball, including the gravitational energy, is still given by the formula 



Etat = / cPxp, (11) 



whereas the energy of the scalar configuration without the gravitational contribution, E s , can be 
defined as 

E s = [ d 3 x VABp. (12) 



Of course, E tot is the true energy of the configuration, but E s is comparable to the corresponding 
flat space configuration, as will be shown. In the presence of gravity also the formula for the 
charge of the Q-ball is modified and is now given by 



Q = -iuj J d 6 xVAB (0V - <f>^) , (13) 

which differs from the flat space formula by the measure factor y/AB. 

From the set of equations (|9|) one can see that the Ansatz (0) for the scalar field is still applica- 
ble. The right hand side of the equation for Rt r then vanishes, implying that A is independent on 
time. Then, as usual, we can without a loss of generality assume that also B is time independent. 
The new set of equations now reads as 
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-AmG{p + 3P)B, 

-4ttG{p-P)A, (14) 



■ 1 + i-hIt-%i = - 4 ^- f >' 2 



and 



The w-dependent potential is defined as 




AV* u (B,<p) = Q. (15) 



V u (B,tp) = -V(tp i )--—<p* (16) 

and denotes a derivate of the potential with respect the field ip. The charge of a Q-ball is now 
given by 

Q = 8ttlo J™ dr r 2 ^A(r)B(r) ip 2 . (17) 

Like in the flat case, a necessary but not a sufficient condition for the existence of Q-balls is 
that the w-dependent potential has a non-zero minimum ip + with V w (B,ip + ) < Vuj(B,0) Q. As 
usual, this restricts the values of u, but it also constrains the function B, because V u depends 
explicitly on it. The allowed values for uj are hence restricted to be larger than a critical, non- 
negative value of u c , i.e. uj > u> c . The critical value, oj c , is determined by the set of equations 

V Uc {l,ip c ) = 0, 

Kfrfc) = o. (18) 

The equations for A and B ( |14] ) can be easily solved in terms of p and P: 

~Y = l-8irGr 2 p, (19) 



1 From now on we suppose, that 14,(0) = 0. We can also choose, without a loss of generality, the central values of A 
and B to be ^4(0) = B(Q) = 1. In addition, for finite energy field configurations A(oo) = 1. 



& + ^ = 8irGAr(p + P). (20) 

These equations can not, however, be integrated straightforwardly because both p and P are 
dependent on A, B and the field p. Therefore our aim in the following sections is twofold: In the 
next section we study the equations and their solutions by utilising analytical approximations 
in order to gain an understanding of their fundamental properties. To verify these results, Eqs. 
(|l~4] ) are studied numerically in Section 



3 The thin-wall limit and gravity 

If the w-dependent potential V w is such that it fulfils the conventional thin-wall requirements, 
i.e. the height of the barrier is much larger than the difference between minima, the field forms 
a bubble of radius R. Inside the bubble the value of the field is essentially constant and outside 
virtually zero with only a thin transition region separating them. The metric components A and 
B are, however, not constants, but monotonic functions inside the bubble. Thus, one has to 



determine their behaviour directly from Eqs. (1£) and (|20|). To determine the properties of such 



a configuration, it is useful to consider the effective action 

= 4ir J™ dr r 2 ^A(r)B(r) (j^V \rf + V u (B(t), <p(r))j , (21) 



from which the equation of motion for tp, Eq. (|15|), emerges. Applying the usual thin- wall 
approach to the action (|2l|), leads to some complications due to the presence of gravity. From 
Eq. ( |i"8| ) we see that two separate sub-cases need to considered, depending on whether the 
potential V is monotonic or it has a non-zero local minimum. Defining <p + by 

V^,<P+) = 0, Vl(l,ip+)<0, (22) 

i.e. <p+ determines the location of barrier maximum, we can now make difference between type 
I and type II thin-wall Q-balls. A thin-wall solution is classified as type I if V(j c (l,(p+) — tjJ 2 p 2 
and type II if V UJc (l, tp + ) 3> w^ip 2 ,, where tp c is determined by Eqs. (|ilf). Type I thin- wall hence 
corresponds to a monotonically increasing potential V whereas type II thin-wall solutions appear 
when the potential has a non-zero local minimum such that it is nearly degenerate with minimum 
at ip = 0. 

In both cases we assume that inside the bubble tp is essentially constant. Its value tp u is 
determined by equation 

Vt,(B(R),<p u ) = 0. (23) 

Note that even though B(r) is varying inside the bubble, the change is relatively small which 
justifies our approximation, B(r) m B(R). 

The value of inside the bubble is hence given by 

Su,,in = 4vr I* dr r 2 VAB V u (B(r), tp(r)). (24) 



o 

Note, that now uj > lo c and thus V w {B{f), tp(r)) < and E^in is negative. Note also that due 
to Eq. (p3|), tpu depends on the bubble radius R. This negative contribution to E w should be 
balanced by the surface term. If A(R) and B{R) differ only little from unity, the wall contribution 
can be approximated by 

S^wall - ^ctR 2 , (25) 

where the surface tension is determined by varying the action with the degenerated potential V w . 
In practice this means that, as a first approximation, we set iv = uj c , B = A = 1, and we finally 
end up with the usual formula for the surface tension: 



yj2V Wo {l,<p) dp. (26) 



Therefore the effective action, reads out as 



£>ui — ^ui,in ^>uj,waU- (27) 

The difference between Eq. ( p7[ ) and the usual thin-wall case lies in the explicit dependence of 
Eq. (|27| ) on the metric components A and B. The next step is then to solve A and B from Eqs. 
( |I^ ) and ([H])- Note, that in this approximation the energy density 

P=^M 2 + ^'| 2 + tV) (28) 

is virtually constant inside the bubble, where one can also neglect the derivative term. 

From this point on the two bubble types corresponding to the different types of thin-wall 
potentials differ from each other and are discussed separately. 



3.1 Type I thin- wall solution 

For the type I solution the terms ^-|v| 2 an d V{(p 2 ) are essentially equal. The relative variation 
of p is also small and hence it is essentially constant. We can thus integrate A out from Eq. ( |T9| ) 
and obtain 

A(r) = i _ ^Gpr* = l - m-Gv^r 2 ' ^ 

Next we note that within the same approximation p + P = ^-\<p\ 2 + ^Iv 9 '! 2 — 2 ^g _ M 2 an d Eq. 
( |20D can be rewritten as 

A' 

B' + —B = 16TrGu 2 ip 2 Ar, (30) 
which can be solved in the the leading order of G, 

B{r) = l + ^GV^l)r 2 . (31) 

We are now able to write down the equations determining the gravitationally coupled thin- 
wall solution for the Q-balls. We use the solution given above for B(r), solve the relation between 



R and ip^ using (|23|) and calculate R from ^rg- = 0. So, in addition to Eqs. ( [29|) and (|3l[), we 



dR 

have equations 



V^(B(R), <p u ) = 0, AttR 2 sJA{R)B{R) V u (B(R), + SnuR = 0, (32) 

which determine the value of uj as well as R and the corresponding value of the field, <*p u - Using 
these values for u, R and tp u we may also calculate the corresponding charge Q from Eq. (ji~7|) . 
If A and B do not differ too much from unity we simply obtain 

Q = y W ^i? 3 , (33) 

where we have neglected the contribution from the wall. 

The new features of the gravitationally coupled Q-balls arise due to the non-trivial dependence 
of the energy inside the bubble to the bubble radius R. Suppose that uj is close enough ui c . The 
scale factor B{R) is a monotonically increasing function, and it may happen that at some point 



the effective w-parameter ui{R) = uj/y / B(R) becomes too small and the Eqs. (32) cannot be 
satisfied. Because R and (p^ depend implicitly on uj, we are now able to define another critical 
value of uj, ujq, as the smallest value where a solution of Eqs. ( |32"| ) can be found. Thus when 
gravity effects are taken into account, Q-balls seem to exist only if u > log which, on the other 
hand, implies maximal a radius R ma x and maximal charge Qmax = y w G^ G ^mtti- The only 
other possibility is, that the singular point of A(R) is hit first. That, however, would mean that 



such a configuration would collapse and form a black hole and hence such Q-balls are not possible 
either. Clearly, the critical values u> m i n , Rmax and Q m ax depend explicitly on the shape of the 
potential and to proceed further one has to define it more precisely. 

To see explicitly how the Q-ball charge is constrained, we study a simple example potential. 
The model potential, V, is constant, m 4 , for tp < p>x, expect near the origin where it tends to 
zero as p 2 . The exact behaviour near the origin is not important in our analysis. For values 
p > pi the potential increases rapidly, e.g. V ~ p n where the power n is sufficiently large. 

Near wg, the minimum of V u is essentially <p\. The metric components are approximately 
given by 

A(R) = l+ l -^GV{p\)R 2 
B{R) = l + ^GV{pl)R 2 

where we have dropped terms higher or equal to 0(G 2 ). Within the same approximation, we get 
from Eq. (]3^ ) an equation for R: 

iTxGm ^ 3 ( 3m 4 _ w 2^2) + i#( m 4 _ Jl^ + 2a = 0. (34) 



Only one of the roots of the cubic equation (34) is physical and is increasing with decreasing lo 



until the critical point, uj = u>g, where the discriminant of the equation vanishes. The maximal 
value of R is ^ ^ 

Rmax = GvrGm^m 4 -^?)) ' (35) 
The corresponding charge and energy of the configuration are 



Qmax ~ 1(P-) 2 ^- (36) 

m m z G 
pi 1 

Emax ~ 2 — , (37) 

m m(j 

where we have assumed that u>qP 2 ~ m 4 and used the approximation a ~ m 2 p\. The energy 
to charge ratio of the maximal Q-ball is E max /mQ max = m/p^. In order to see that our 
approximation is valid, we note that Gm 4 R 2 nax ~ (pi/Mp\) 2 / 3 1. 

3.2 Type II thin-wall solution 

For the type II solution the fine tuning V U (B, 4>) ~ is not needed in order to reach the thin wall 
limit because a large barrier exists even without the cu 2 (fi 2 -term. Now, indeed, uP'p 2 ^ 3> V(p 2 } ) 
and we write an effective energy density as p = (1 + n)ui 2 p> 2 , where k is assumed to be a constant. 
Instead of equation ( |29| ) we now write 

A ^ = i _ ^Gpr 2 = l-^-G{l + K)uj 2 p 2 r 2 (3§) 
and from Eq. (|30| ) we obtain by integrating 

B(r) = 1 + ^G(2 - K)u?p 2 r 2 . (39) 



1 



The difference between type I and II solutions is that in the type I solution, A increases near 
the origin twice as fast as B whereas for the type II solution the roles are reversed: B increases 
roughly twice as fast as A whenever k<1. 



Due to exactly same reason as before, a minimal co appears together with a maximal bubble 
radius, charge and energy in this case as well. Generally, they can again be found solving Eqs. 
(|32|), ( |38D and (^9[) and searching for the smallest possible oj. However, the procedure is strongly 
dependent on the explicit form of the potential V as it was in the case of type I solutions. 

Therefore, as an example, we again proceed by studying a simple potential 

V(ip 2 ) = \cp 2 (((/? - Mf + eM 2 ) , (40) 

where M is a mass scale and e <C 1 is a small positive parameter lifting the non-zero minimum 
up from V = 0. Note, that the quadratic term near the origin reads as A(l + e)M 2 ip 2 , so that 
the mass of a single quantum is given by m 2 ~ XM 2 . Calculating the surface energy, we find 
that a = \/AM 3 /6. The set of equations (|32^) reads now as: 

\{^-M){2^-M)-^ )+ u 2 c = 0, (41) 

UJ 2 , n\ 2(7 



A(R)B(R)(\vl(^-M) 2 + (u J 2 c -—)vi)+— = 0, (42) 
where uj 2 = XeM 2 . 

Evidently uiq is close to uj c and (p w is close to M, which allows us to write ip w = M and 
uj = uj c wherever no explicit difference appears in the formulae. Expanding in the first power of 
G we arrive at an equation for the bubble radius, 

3 Sj^/u; 2 - 1) 3a 

87tGlu 2 M 2 {2 - k) 2ttG^M 4 (2 - «) 1 ; 

Following the same arguments as in the type I case, we can solve for R max and w^: 

R ~ ( ^ Mg M V3 (44) 
Kmax ~ W{2-n)Mui) [ ' 

u; 2 G = i o 2 (l + 2((2- K )n-^) 1 / 3 ). (45) 

The approximations we made in the process of deriving the above equations constrain the ac- 
ceptable range of parameter values, M <C y/Mpyuj c . If this condition is not satisfied our results 
are not valid and the equations need to be solved numerically. 
The maximal charge and energy of a Q-ball in this scenario are 



2VX MM 2 , 

3(2 - K) UJ* 

2VA MM 2 , 



The energy to charge ratio, E max /mQ max = ui c /m = y/e/ <C 1, indicates that the maximal Q-ball 
is energetically stable. 



4 Numerical Results 

Next we turn to look at some numerical results obtained by computer simulations. We concentrate 
ourselves to the two different potentials, one with a flat plateau and one with a nearly degenerate 
minimum. The flat potential allows type I thin wall Q-ball solutions near the critical case. Such 
potentials are hence called type I potentials. The other type of potential is the one given by 



Eq. (|40|), leading to the type II thin- walls solutions and are therefore called type II potentials. 
Note, that both of these potential types are represented, e.g., in the MSSM, where the potential 
arising from gauge mediated supersymmetry breaking is clearly a type I potential whereas gravity 
mediated supersymmetry breaking typically leads to a type II potential. 

From the simulations we have learned, however, that the thin wall limit is exceptionally 
difficult to study by simulation due to an extreme fine tuning required from the initial value of 
the field inside the bubble. This fine tuning is much more severe in the gravitationally coupled 
case, because here one is not attempting to find just the field value, but also a corresponding 
value of u). Therefore the simulations are concentrated to the thick- wall end of the parameter 
space, i.e. to larger oj. At best we are able to to see some hints of the beginning thin-wall 
behaviour. Also, in order to clearly demonstrate the observed behaviour of the solutions, we use 
large mass scales in the potential, since otherwise the true thin-wall solution is needed. 



4.1 Type I potential 

We begin with a potential 

{m 2 ip 2 , \ip\ < ip 

m 2 <p 2 , ip Q < \(p\ < if! (48) 

™Vo + e (^ 4 - <Pi), <Pi < <P 

where m 2 determines the curvature at the origin, e is a small positive number and < (po <C (pi- 
To ensure that the minimum is in the region <p > ip\ and to is larger that the critical value, 
we have to require m 2 > lo 2 > 2etp 2 and 4e(m 2 </?Q — eyf) < w 4 , respectively. This potential 
resembles the flat one discussed in the context of thin- wall configurations. As an example we 
choose m = ipo = 10 _3 Mpi, ipf = ilO _8 Mp[, and e = 10 -4 . For these values, the critical oo 
equals u> c = y/2 X 10 _1 m and the critical field value is just ip\. The configuration is studied for 
several values of u in order to demonstrate the dependence of the solution on u. Hence, we take 
uj/m = 0.5, 0.4, 0.3, 0.2 and 0.142. In Fig. [j]a we have plotted the field configuration of p and in 
Fig. the ratio (B(r) — l)/(A(r) — 1) for these values. Furthermore, in Fig. || we have plotted 
the corresponding energy of the free quanta mQ, total energy E tot and scalar field energy E s as a 
function of r for the largest Q-ball corresponding oj/m = 0.142. Note that in Fig. [2], the energy 
of the free quanta is scaled by a factor of ten, i.e. Etot, E s <C mQ. 
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Figure 1: (a) The field profile as function of radial coordinate, (b) The ratio of metric compo- 
nents (B — 1)/(A — 1). In the both figures curves represent values ujjm = 0.5 (the 
shortest curve), 0.4, 0.3, 0.2 and 0.142 (the longest curve). 
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Figure 2: The scaled energy of free quanta 0.1 x mQ/m (solid line), the total energy E^at/rn 
(dashed line) and the scalar field energy E s /m (dotted line) inside radius rm. 



Evidently, for large u, the profile of cp is thick- walled. However, from the figure it is clear that 
when u decreases the configuration tends towards the thin-wall case. This is also reflected by the 
fact that the ratio (B(r) — l)/(A{r) — 1) decreases along with decreasing u. Indeed, according to 
our analysis of the type I thin- wall configurations, inside the bubble the ratio should approach i. 
As can bee seen from Fig. ^ the energy difference between E to t and E s is already for uo/m = 0.142 
about 7% and evidently it is increasing with increasing Q-ball size. It should be noted, that the 
critical to which is now oj c /m ~ 0.1414, is already quite close to the smallest value used in plotting 
the figure. In the large Q-ball limit, however, the configuration is extremely sensitive to the value 
of OJ. 



4.2 Type II potential 
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Figure 3: (a) The field profile as function of radial coordinate, (b) The ratio of metric compo- 
nents (B — 1)/(A — 1). In the both figures curves represent values tv/m = 0.5, 0.4, 0.3 
and 0.261. 



For numerical studies of potential (^0| ) we choose A = 1, M = 0.0316Mpi and e = 10 3 . 
In Fig. ^a we have plotted the field configuration corresponding to several different values of 
u/m = 0.5, 0.4, 0.3, 0.261. The ratio (B(r) - l)/(A(r) - 1), which inside the bubble equals 
(2 — k)/(1 + k) , is presented in Fig. ||b. Moreover, we now also present in Fig. [| the energies -Etot 
and E s as function of the charge Q. For these values of m, A and e the critical uj is simply given by 



= 10 _3 m 2 . Note, that this time the smallest u> we are able to calculate is oj = 0.261m which 
is significantly bigger that the critical cj — 0.0316m. This reflects the fact, that the thin- wall 
type behaviour begins at much larger ratio uj/uj c for a type II potential that for a type I potential. 
From Fig. ^b one can also read that the effective parameter k inside the bubble for small u> is 
about 0.07. The difference between total and scalar field energies, for the largest plotted balls 
Q = 4.8 x 10 3 , is about 13%. 
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Figure 4: The total energy 2? tot /m (solid line) and the scalar field energy E s /m (dotted line) 
as function of the charge Q. The smallest Q corresponds to uj/m = 0.995 whereas 
the largest one corresponds to uo/m = 0.261. The straight line is the E = mQ curve. 



5 Discussion and Conclusions 

In this paper we have considered the effects of gravity on Q-ball type scalar field configurations. 
Including the effects of gravity leaves the flat space Q-ball Ansatz valid; one still has a spherically 
symmetric solution whose complex phase rotates uniformly with time. The geometry of space- 
time is modified from the flat case, which gives a contribution to the charge and energy of the 
configuration. 

From the analytical and numerical work it is evident that the effects of gravity can give a 
large contribution to the Q-ball energy. However, in order for the contribution to be significant, 
either the free quanta of the Q-ball scalar field need to be extremely heavy or the charge of the 
Q-ball needs to be very large. The gravitational contribution is strongly dependent on the details 
of the potential and needs, in general, to be computed numerically. 

Interesting points are raised by the analytical analysis. The analytical work is concentrated 
on the thin- wall limit, in order to gain an understanding of the properties of large Q-balls for 
which gravitational effects are more likely to be significant. In the thin-wall limit, analytical 
considerations reveal that the maximal size of Q-balls is ultimately limited by gravity. The 
non-trivial geometry inside the Q-ball limits the size of Q-balls for any potential that allows for 
Q-ball solutions. The physical reasoning behind this is simple. In equilibrium, the rotation of 
the complex phase generates an outward pressure which prevents the Q-ball from collapsing due 
to gravitational forces. In the ordinary thin-wall limit, the Q-ball frequency tends a critical value 
u c , while the charge the Q-ball becomes larger and larger. If gravity is included in the scenario, 
it is evident that at some point the Q-ball becomes too large to be balanced by the internal 
pressure. Hence, there exists another critical value, log > oj c , corresponding to the maximally 
large Q-ball. 

The limit on the maximal Q-ball size due to the gravitational effects is an interesting point 
raised by this analysis. The two example potentials discussed approximately correspond to the 



potentials associated with gauge (type I) and gravity (type II) mediated supersymmetry breaking. 
Using typical values for the variables, m ~ 100 GeV, ipi ~ M ~ 10 15 GeV, and approximating 
lu c ~ m, estimates for the maximum charge, energy and radius are: Qmax ~ 10 60 , E mSLX ~ 
10 49 GeV, R max ~ 10 13 GeV" 1 - 0.1cm (type I) and Q max ~ 10 34 , E max ~ 10 36 GeV, R max ~ 
1 GeV" 1 (type II). In both cases, the Schwarzschild radius of a corresponding black hole is still 
much smaller than that of the Q-ball. Gravity will clearly not have a very significant role to play 
in a cosmological setting for Q-balls in the gauge mediated scenario. In the gravity mediated 
scenario the limit is much lower and may be significant. It is interesting to note that, in order 
for the primordial Q-balls to act as a self- interacting dark matter candidate [16|, one typically 
has to consider Q-balls with very large charges fll7fl . 

Since the gravitational contribution to the Q-ball energy is negative, it is possible that gravity 
can render an otherwise energetically unstable Q-ball stable. This behaviour was also confirmed 
by the numerical work, where is was observed that the negative energy of the gravitational field 
could stabilise a large enough Q-ball. One can then envision a cosmological scenario where Q- 
balls are unstable, and hence do not exist, unless their charge is large enough for the gravitational 
effects to stabilise them. A primordial Q-ball distribution formed at a large energy scale would 
then leave behind only stable, super heavy relic Q-balls. 
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